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ABSTRACT 

We show that a star orbiting close enough to an adiabatically grown supermassive black hole can capture a 
large number of weakly interacting massive particles (WIMPs) during its lifetime. WIMP annihilation energy 
release in low- to medium-mass stars is comparable with or even exceeds the luminosity of such stars due to 
thermonuclear burning. The excessive energy release in the stellar core may result in an evolution scenario 
different from what is expected for a regular star. The model thus predicts the existence of unusual stars within 
the central parsec of galactic nuclei. If found, such stars would provide evidence for the existence of particle 
dark matter. White dwarfs seem to be the most promising candidates to look for. The signature of a white 
dwarf burning WIMPs would be a very hot star with mass and radius characteristic for a white dwarf, but with 
luminosity exceeding the typical luminosity of a white dwarf by orders of magnitude <50L Q . A white dwarf 
with a highly eccentric orbit around the central black hole may exhibit variations in brightness correlated with 
the orbital phase. 

Subject headings: black hole physics — elementary particles — radiation mechanisms: non-thermal — stars: 
general — stars: evolution — dark matter 



1. INTRODUCTION 

The nature of the non-baryonic dark matter, which dom- 
inates the visible matter by about 4:1, is perhaps the most 
interesting experimental challenge for contemporary particle 
astrophysics. A hint for a solution has been found in parti- 
cle physics where the WIMPs arise naturall y in supersym- 
metri c extensions of the Standard Model (e.g jHaber & Kariel 
fl985h . among other possibilities. The WIMP is typically de- 
fined as a stable, electrically neutral, massive particle. Assum- 
ing that non-baryonic dark matter is dominated by WIMPs, 
the pair annihilation cro ss-section is restricted to provide 
the observed relic dens ity (lJungman. Kam ionkowski & Griestl 
1996; Bergstrom 20(3). A pair of WIMPs can annihilate pro- 
ducing ordinary particles and 7-rays. 

WIMPs are expected to form high density clumps ac- 
cording to N-body simulat ions of test particles with only 
gravitational inter actions dNavarro. Frenk. & White! 119971 : 
iMoore et all Il999t) . The highest density "free space" dark 
matter regions occur for dark matter particles captured within 
the gravitational p otential of adiabatically grown supermas- 
sive b lack holes dGondolo & SilkJ 119991: iBertone & Merrittl 
2005). Higher dark matter densities are possible for dark mat- 
ter particles captured inside of stars or planets. Any star close 
enough to an adiabatically grown supermassive black hole can 
capture a large number of WIMPs during a short period of 
time. 

Such an idea has been fir st proposed byl Salati & Silk ( 1989) 
and further developed by Bouque r& Salatil d 1989b who ap- 
plied it to main-sequence stars. WIMP annihilation in stars 
may lead to considerable energy release in the stellar cores 
thus affecting the evolution and appearance of such stars. The 
model led to the conclusion of suppression of stellar core con- 
vection, thus predicting a concentration of stars in the Galactic 

1 Also Kavli Institute for Particle Astrophysics and Cosmology, Stanford 
University, Stanford, CA 94309 



center masquerading as cold red giants. 

We perform order-of-magnitude estimates of WIMP cap- 
ture rates for stars of various masses and evolution stages 
located in high density dark matter regions. We use current 
limits on WIMP-nucleus interaction and WIMP annihilation 
cross sections, as well as recent estimates of WIMP energy 
density near an adiabatically grown supermassive black hole. 
We argue that white dwarfs, fully burned stars without their 
own energy supply, are the most promising candidates to look 
for. 

2. WIMP ACCUMULATION IN STARS 

In a steady state t he WIMP capture rate C is balanced by 
the annihilation rate dGriest & Seckellll987t) 

C = AN 2 X (1) 

where N x is the total number of WIMPs in the star, 

A = (cy a v)/V cS , (2) 

V cif = Tr 3/2 r 3 x , (3) 

(cr a v) is the velocity averaged WIMP pair annihilation cross- 
section, and the effective volume V c s is determined by match- 
ing the core temperature T c with the gravitational potential 
energy at the core radius r x (assuming thermal equilibrium). 
The total number of captured WIMPs is 

N x = Cr e gtaith(T*/T e g), (4) 

where r* is the star age, and the equilibrium time scale is 
given by 

T eq = {CA)- x ' 2 . (5) 

The numbe r density distribution o f WIMPs can be es- 
timated as dPress & Soergell Il985t iGriest & Seckel 119871: 
iBottino et al.ll2002l) : 



U x( r ) = K eX P(- r2 Ay)> 



(6) 
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FIG. 1 . — Mass-radius-capture rate diagram for Hydrogen stars. The dotted 
line series are calculated using eq. (9) and the numbers on the plot show the 
logarithm of WIMP capture rate log 10 C. The strong dash ed line labeled 
"MS" shows characteristics of the main sequence stars ( Massey & Mever 
20011) . and the diamond symbol indicates a sun-like star with M = Mq, 
R = Rq . The open circles on the main sequence show the stars of particular 
spectral type and the numbers show their luminosity log 10 (L t / Lq). The 
solid line divides the diagram into "allowed" and "forbidden" parts, see text 

w"h m fe e n^ tal = N x /V e s is the central WIMP number density. 
In thermal equilibrium, the effective radius r x = tt is deter- 
mined by the core temperature T c and density p c 



r T 



( 3T C 



\2 / nGp c m x 



1/2 



(7) 



where c is the speed of light, G is the gravitational constant, 
and m x is the WIMP mass. 

If the WIMP-nucleon scattering cross-section is large 
enough, then every WIMP crossing the star is captured and C 
saturates to a maximal value proportional to the star's cross- 
sectional area 7ri? 2 , the WIMP velocity dispersion v, and the 
ambient WIMP energy density p x : 



C = 



3?r 



1/2 



Px v 
m x 



^ 2v 2 



(8) 



where £ is taken to be 1 .77, v esc = \/2GM*/ R* is the escape 
velocity at the star's surface, and Af* is the star mass. 

On the other hand, limits from direct detection of dark 
matter on the WIMP-nucleon cross-section imply that only 
a fraction of the WIMPs crossing the star will scatter and 
be captured. In this case the star's cross-sectional area is re- 
placed by the product of the WIMP-nucleon scattering cross- 
section cr„, the number of nuclei, and a spin-independent 
coherent factor for nuclei heavier than Hydrogen, that is 
A\ for a nucleus with the atomic mass A n . The limits on 
spin-independent WIMP-nucleon scattering are more strin- 
gent than for spin-dependent; from direct detection experi- 
ments the spin-independent WIMP- nucleon scatte ring cross- 
section cr is at most 1(T 43 cm 2 dCDMSl 12006b . The cur- 
rent limit on spin-dependent Hydrogen scattering comes from 
Super- Kamiokande, which finds that a s is at most 10~ 38 cm 2 
( Super- Kamiokande 2004). In the following formulas we will 
use a generic cross-section <j n to be substituted by ao or a s 
for spin-independent and spin-dependent interactions, corre- 
spondingly. If the star is entirely composed of nuclei with the 
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FIG. 2. — Mass-radius-capture rate diagram for Helium stars. Lines are 

a¥omfc ma'sl 111 ^ (A n — 1 for pure Hydrogen), we have for 

m x ^ A n m p : 

1/2 
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(9) 



where mp is the proton mass. 

Eq. ([Sf gives an absolute upper limit on the capture rate. 
The allowable range of masses and radii at a given WIMP- 
nucleon scattering cross section can be obtaned from eqs. ^ 
and ©: 



R* > Al!\ 



imip 



(10) 



Smaller radii are forbidden as the capture rate will exceed the 
upper limit given by eq. ©. 

A star spends about 90% of i ts life burning hydro gen and 
most of the rest burning helium dWooslev et al.ll2002l and ref- 
erences therein); the carbon, neon, oxygen, and silicon burn- 
ing lasts about 1000 yr altogether. The star core temperature, 
matter density, the radius, and the mass, all change during the 
evolution cycle. The mass-loss by a star during its lifetime im- 
plies that stars with initial mass > 35Af and the solar metal- 
licity to end their life as hydrogen-free objects of roughly 
5M ; early massive stars with essentially no metallicity will 
evolve differently and may retain a massive Hydrogen enve- 
lope even at later burning stages. 

Figures QJ[3] show the mass-radius-capture rate diagrams 
for m x = 100 GeV, (a a v) = 3 x 10~ 26 cm 3 s~\ and 
WIMP velocity dispersion v = 270 km s _1 . The value 
p x — 10 10 GeV cm -3 corresponds to the maximal en- 
ergy density allowed by the age of the supermassive black 
hole ~10 Gyr, and our selected values of {og v) and m x 
(lGondolo&Silklll999l: iBertone & Merrittll2005l) . The spin- 
dependent and spin-independent WIMP-nucleon scattering 
cross sections are taken as a s = 10~ 38 cm 2 and errj = 10~ 43 
cm 2 , correspondingly. The series of dotted lines are calculated 
using eq. (O, and the numbers on the plots show the logarithm 
of WIMP capture rate log 10 C. The strong dashed line labeled 
"MS " in Figured] shows char acteristics of the main sequence 
stars (Massey & Meyer 2001), and the diamond symbol indi- 
cates sun-like star with M = Mq, R = R Q . The open cir- 
cles on the main sequence show the stars of particular spectral 
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TABLE 1 

WIMP ACCUMULATION DURING DIFFERENT BURNING STAGES 
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NOTE. — Af init is the initial hydrogen star mass. 
d Stellar model with 10 -4 solar metallicity and assuming a massive Hydrogen envelope (<r n — a s ). 
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FIG. 3. — Mass-radius-capture rate diagram for Carbon stars. The rectan- 
gular area outlined by the dashed line and labeled "WD" shows parameter 
range typical for white dwarfs. Lines are coded as in FigurefT] 

type and the numbers show their luminosity log 10 (L»/L Q ). 
The rectangular area outlined by the dashed line and labeled 
"WD" in Figure[3]shows the parameter range typical for white 
dwarfs. The solid line (eq. fTOl ) divides the diagram into "al- 
lowed" and "forbidden" parts. For any particular M» and 
the WIMP capture rate can be found from the grid provided 
by the dotted lines. 

Table Q] shows the WIMP capture rate, WIMP accumu- 
lation, and other parameters for different star masses (1 — 
75M Q ) and at different evolution stages. Spin-dependent 
WIMP accumulation is calculated for the Hydrogen burn- 
ing stage and spin-independent WIMP accumulation is cal- 



culated for the Helium and Carbon burning stages. The star 
masses, radii, c ore te mperatures and densities are taken from 
IWoosley et alJ {2002) and correspond to the particular burn- 
ing stage. The r* column shows the duration of the corre- 
sponding burning cycle. The uppermost row in the Carbon 
burning stage table shows the parameters typical for a white 
dwarf where the core temperature is taken ad hoc. This does 
not affect the WIMP capture rate, however. Smaller annihila- 
tion cross section ((J a v) < 3 x 10~ 26 cm 3 s _1 will lead to 
a larger number of captured WIMPs, although this would not 
change the burning rate. Early massive low-metallicity stars 
can potentially accumulate the most numbers of WIMPs due 
to the presence of a large Hydrogen envelope. Such stars could 
capture WIMPs at greater rates due to the potentially larger 
spin-dependent WIMP-nucleon scattering cross-section. The 
WIMP number density may f urther increase du ring the sud- 
den collapse of a massive star (lUllio et al.ll2.Q01h . 

In steady state the annihilation rate is equal to the capture 
rate. The energy release due to the WIMP annihilation in the 
star core is 



L x ~ 1.6 X 10 34 C 3E 



erg s 



(11) 



.100 GeV 

where C35 is the capture rate in units 10 35 s . 

3. DISCUSSION AND CONCLUSION 

Where does the energy released during the WIMP annihi- 
lation go? Table Q] shows that the typical radius of the ther- 
mal distribution of WIMPs in the star core rth *C R*, there- 
fore, the products of their annihilation can not propagate to 
the star surface and are eventually converted into the thermal 
energy and neutrino emission. This may be not true for the 
white dwarfs whose radii are comparable to r t h, and where 
the effect of gravitational "focusing" is important (eqs. (8), 
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|9l , second term in square brackets). While most of the en- 
ergy is converted into heat, the surface of white dwarfs may 
emit products of WIMP annihilation, particles (e , p, p, v, v) 
and 7-rays, and this might be visible with a 7-ray telescope. 

The WIMP capture rate generally increases with the star 
mass. However, the relative importance of WIMP annihila- 
tion diminishes in massive stars since their own luminosity 
increases as oc 5 . Therefore, the interiors of massive stars 
will remain largely unaffected by WIMPs. 

The WIMP burning affects mostly stars with M < M 
with low to moderate luminosities and a large rate of WIMP 
capture. The typical luminosity of a sun-like star burning Hy- 
drogen is L = 3.8 x 10 33 erg s _1 , while WIMP annihilation 
in its core can provide up to L x ~ 4.6 x 10 34 erg s _1 , i.e. 
10 times more. Such an additional energy source may exceed 
the star's own resources from nuclear burning. The large con- 
centration of WIMPs may a lso lead to more effectiv e energy 
transport in the stellar core (Bouquet & Salati 1989) thus de- 
creasing the core temperature and limiting the replenishment 
of the burning core with fresh nuclear fuel. Combined, these 
two effects could change the star's evolution scenario, life ex- 
pectancy, and appearance. 

A white dwarf, a star without its own energy supply consist- 
ing of Carbon and Oxygen, may emit L x ~ 2 x 10 33 — 2 x 10 35 
erg s _1 , i.e. 0.4-50 times luminosity of the sun, burning 
WIMPs only and this energy source will last forever! To reach 
such a luminosity, the surface temperature of the white dwarf 
should be in the range of T ~ (80 - 150) x 10 3 K assum- 
ing R = O.Oli?©. In comparison, the luminosity of a regular 
white dwarf does not exceed ~0.01L Q . The interiors of white 
dwarfs are almost isothermal, and t he energy trans port is dom- 
inated by degenerate electrons (see Han senll2004 for a recent 
review); therefore, the large number of captured WIMPs and 
their annihilation in the core would not change the structure 
of white dwarfs. Interestingly, since L x oc p x , a population 
of such white dwarfs located at different distances from the 
central black hole will exhibit a luminosity correlated with 
the radial WIMP density distribution. A signature of a "dark 
matter burner" would be a very hot star with a mass and ra- 
dius characteristic for a white dwarf, but with a luminosity 
exceeding the typical luminosity for a white dwarf by or- 
ders of magnitude <50L Q . Note that an independent deter- 
mination of the M*/li* ratio is possible using the gravita- 
tional redshift which has to b e equivalent to vr ~ 50 km s _1 
dGreenstein & Trimbldfl967t) . 



Furthermore, a low-mass star with a highly eccentric orbit 
around the central black hole may exhibit variations in bright- 
ness correlated with the orbital phase. The brightness should 
increase as the star approaches the periastron and some time 
after that. To have this working, the orbital period should es- 
sentially exceed r eq . Carbon burning stars have r eq ~ 10 yr. It 
is even shorter in case of a white dwarf r eq ~ 0.5 yr. If a white 
dwarf appears in a high-density WIMP region, the WIMP 
density in its material quickly reaches equilibrium. On the 
other hand, the small radii of white dwarfs i?* ~ lQrth im- 
ply that the captured WIMP density close to the surface may 
be large enough thus the surrounding WIMP density change 
will result in variations in their brightness. This makes white 
dwarfs ideal objects to test the WIMP density in the environ- 
ment in which they are orbiting. 

Advances in near-IR instrumentation have made possible 
the observation of stars in the inner parsec of the Galaxy 
dGenzel et all l2000t iGhez et alj l200l 12001 . The observed 
absorption line widths imply high temperatures and lead to 
a "paradox of youth:" apparently young stars in the region 
whose current conditions seem to be inhospitable to star for- 
mation. One of the possibilities is that they are old stars mas- 
querading as youths. If an independent determination of their 
mass reveals that some of them are low-mass stars, then they 
would become candidates for the dark matter burners. 

The model predicts the existence of unusual stars and con- 
sequently unusual supernovae in the central parts of the galax- 
ies. These are the stars burning WIMPs at a high rate whose 
luminosities essentially exceed the typical luminosity of a star 
of given mass. If found, such stars would be interesting probes 
of particle dark matter near supermassive black holes. Their 
luminosity, or rather their excess luminosity, attributed to the 
WIMP burning can be used to derive the WIMP matter den- 
sity at their location. On the other hand, a lack of such un- 
usual stars may provide constraints on WIMP density, WIMP- 
nucleus scattering and pair annihilation cross-sections. 
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